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Riassunto. Viene definito un operatore LG, costituito da differenze finite, calco- 
late su un insieme discreto G C R, che approssima l’operatore di Kolmogorov 


L=02+r0,-0d:. 


La viene costruito sostituendo un rapporto incrementale secondo alla derivata e, 
ed un rapporto incrementale lungo il vettore v = (0,x,—1) alla derivata prima 
x, — d;. Viene poi fornita una condizione di stabilità per le soluzioni di Lew = 0. 

Successivamente si studia un problema di Dirichlet relativo all’equazione di Kol- 
mogorov e si dimostra che, se vale la condizione di stabilità, le soluzioni del problema 
di Dirichlet relativo ad Lg convergono alla soluzione del problema di Dirichlet re- 
lativo ad L. 


Abstract. We define a grid G C R3 and a “discrete” operator Le, which is 
consistent with the Kolmogorov operator 


L=02+r0,-d. 


We obtain such an operator by substituting the derivatives of the operator L with 
finite differences. In particular the first derivative rd, — è is replaced by a finite 
difference along the vector v = (0,x,—1). 

Next we give a stability condition for the solutions of Lew = 0, which is the 
same condition corresponding to the heat equation in two variables: 02 — 0;. Under 
this condition we prove a maximum principle relative to Lc. 

Then we consider a Dirichlet problem for the Kolmogorov equation Lu = 0 and 
define a corresponding “discrete” Dirichlet problem for Lew = 0. Using the above 
maximum principle we show that the solutions of the discrete Dirichlet problem 
converge to the solution of the Dirichlet problem relative to Lu = 0. 
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1. Introduzione e principali risultati. 


I risultati che esporrò in questo seminario sono stati ottenuti in collaborazione 
con C. Mogavero, e costituiscono una parte della sua Tesi di Laurea [Mo]. 

Presenterò un metodo di approssimazione numerica per la soluzione di un pro- 
blema di Dirichlet relativo all’equazione di Kolmogorov in R3 


(1.1) Lu= zz +cu,—w=0 


e all’insieme 
Q= {Gu eR°:u1<2<1,-1<y<10<6<1) 


L'equazione .di Kolmogorov (1.1) è l’esempio più semplice della seguente classe 
di equazioni ultraparaboliche degeneri definite per (2,6) = (21,..., cn, t) e RN+H. 


N N 
(1.2) bp Qi,jUx;,x; + È, di,jLiUx; -—u=0, 


t,j=1 i,j=1 


dove A = (a;;) e B = (b;,;) sono matrici costanti N x N,A è simmetrica e 
semidefinita positiva ed è verificata l’ipotesi: 


[H1]. Tra le righe delle matrici A, AB,.... AB*,... ve ne sono N linearmente in- 
dipendenti. 


Queste equazioni sono dette equazioni di tipo Kolmogorov-Fokker-Planck, e 
godono di proprietà di regolarità analoghe a quelle degli operatori parabolici. 

La proprietà più notevole è senz'altro la seguente: se indichiamo con (2,6) la 
soluzione fondamentale di (1.2) con polo nel punto (e RMH, allora T(z,0) è di 
classe C° in RV+1\{c}. 

Queste proprietà sono strettamente collegate ad una struttura di gruppo di 
traslazioni (non euclidee) e ad un gruppo di dilatazioni, rispetto alle quali gli ope- 
ratori suddetti sono invarianti. Nel caso particolare dell’operatore (1.1), per ogni 
(€É, n.7) € R$, la traslazione Teé,n,r) è definita da 


(1.3) Té,nr)(£,y,t) = (c+É,y+n-té,t+7), 
e, per ogni A > 0, la dilatazione D(A) è 
(1.4) D(A)(z,y,t)= (Az, A5Y, X°t) . 


Per questi risultati e per le motivazioni che hanno portato a studiarli rimandiamo 
senz'altro al lavoro [LP] ed alla bibliografia ivi citata. Ricordiamo solamente che 
equazioni del tipo (1.2) intervengono nella descrizione del moto di particelle brow- 
niane in un fluido. 

Per quanto riguarda il problema di Dirichlet relativo all’equazione di Kolmogorov 
sull’insieme Q, notiamo che non è chiaro a priori su quali punti di 0Q sia naturale 
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assegnare il dato. Per questo motivo abbiamo dapprima considerato il problema di 

Dirichlet “generalizzato” 
Lu=0 in 

(1.5) I Q 


u=g in00Q, 


con g € C(Q), nel contesto della teoria degli spazi armonici astratti. Innanzitutto 
si determina una “soluzione generalizzata nel senso di Perron-Wiener” u € C°°(Q) 
del problema (1.5). Tale funzione risolve in Q l’equazione Lu = 0, ma non è detto 
che assuma il dato < in tutti i punti di 00. 

Il problema di individuare i punti regolari del bordo di un aperto, nei quali il dato 
viene effettivamente assunto, è stato studiato in [Ma] nel caso più generale degli 
operatori (1.2) e per aperti limitati qualunque. Nel caso che stiamo considerando, 
si ottiene una caratterizzazione completa dell’insieme 0,Q dei punti regolari di 00, 
che descriveremo dettagliatamente nel paragrafo 2. 

I risultati sopra richiamati ed un principio del massimo per le soluzioni di Lu = 0 
assicurano l’esistenza e l’unicità di una funzione u € C(Q) N C*(Q) soluzione del 
problema di Dirichlet “classico” 


(1.6) 


Lu=0 inQ 
u=% in0;0Q. 


Non riportiamo qui la prova di quest’ultima affermazione. 
Nel paragrafo 3 costruiremo un operatore Lc, costituito da differenze finite, 
calcolate nei punti dell’insieme 


(1.7) G= {(jAx,kKAy,nAx) e R:j,k,n € z}, 
che sia consistente con l'operatore L, nel senso che 
(1.8) Lgu= Lu+ ((Az)° +A4:)O(Az, Ar). 


Qui O(A:, A:) indica una funzione limitata in un intorno dello zero. Osserviamo 
che per garantire la (1.7) sarà necessario imporre che 


(19) fire 

Potremo inoltre osservare che l’insieme 
H=GNQ, 

su cui studieremo il problema di Dirichlet discreto 


Lew=0 in Ho 
(1.10) 


w=g ind,H, 
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si divide in maniera naturale in due sottoinsiemi: OrH su cui si assegna il dato, e 
Ho in cui l'operatore Lg è ben definito, nel senso che mette in relazione solamente 
i punti di H. 

Nel paragrafo 4 si dimostra un risultato di convergenza delle soluzioni del pro- 
blema (1.10) alla soluzione di (1.5). 

Proveremo innanzitutto che, se supponiamo 


Ai 1 
(1.11) IUNSE, “> 


allora vale un principio di massimo per le soluzioni del problema (1.10) (Teorema 
4.1). 

La (1.11) viene detta condizione di stabilità, in quanto assicura che le soluzioni 
approssimate costituiscono una famiglia limitata (Corollario 4.2). 

Osserviamo che la condizione (1.11) è esattamente la stessa che si richiede per 
l'equazione del calore in due variabili u,, = Ut, ricordiamo però che nel nostro caso 
deve valere anche la (1.8). 

Una griglia G tale che 


(1.12) (Ax, Ayy Ax) = (A, (3) 


ad esempio, verifica le condizioni richieste. Notiamo l’evidente analogia con il 
gruppo di dilatazioni (1.4). 

Infine, utilizzando il principio del massimo e la consistenza dell’operatore LG, 
dimostriamo un risultato di convergenza delle soluzioni approssimate, sotto una 
forte ipotesi di regolarità della soluzione di (1.5) (Teorema 4.2). Tale ipotesi sup- 
plementare può in effetti essere rimossa, ma non riportiamo in questa sede la prova 
completa del risultato. 

Vogliamo infine osservare che il problema (1.10) si presta facilmente ad essere 
tradotto in un programma ricorsivo esplicito da eseguire con un calcolatore. I grafici 
allegati costituiscono il risultato di alcuni esperimenti che abbiamo effettuato su un 
personal computer. 


Ringraziamenti. Desidero ringraziare V. Simoncini per le interessanti discussioni sugli aspetti 
numerici di questo problema e per l’aiuto nella realizzazione degli esperimenti. 


2. Problema di Dirichlet generalizzato. 


Riportiamo ora, senza dimostrazione, alcuni risultati relativi al problema di Diri- 
clet generalizzato (1.5). Introduciamo alcune notazioni: 


Definizione 2.1. Un punto 2 € 0Q si dice regolare per (1.5) se, qualunque sia 
p E C(Q), risulta 


(2.1) lm u(7) = 9(2). 


x€Q 


Chiamiamo bordo regolare di Q l’insieme d-Q dei punti regolari per (1.5). 
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Definizione 2.2. Dato z € 0Q, diremo che il vettore n = (n1,n2,n3) è una 
normale esterna in z se esiste p > 0 tale che 


B(z+pn,p)NQ=0, 


essendo B(x,p)= {ye R3:|c-y|< p}. 


Indichiamo 
er (1,0, 0), v= v(£, Y t) a (0, 2, 1), 


e definiamo, secondo la classificazione introdotta da Fichera in [F], 


x0 = {(7,y,8) € 0Q : (e1,n) = 0 Vn normale esterna in (7,9,t)}, 
Lo = {(,y,9) e DO : (v,n) = 0Vn normale esterna in (2,9,t)}, 

(2.2) Zi = {(,y,) € L° : 3n normale esterna in (7,y,t), (v, n) < 0}, 
Lo = {(2,y,4) € L° : 3n normale esterna in (2,y,t), (v, n) > 0}, 
3 = 00\2°. 


Possiamo verificare direttamente che 


Do ={(0,1,):0<#<1}u{(0,-1,4):0<t<1}, 

Li ={(,y,):-1<2<1-1<y<1}U{(,1):-1<2<0,0<t< 1}u 
{18 :0<2< 1,0<t<1}, 

th = {(c,y,0):-1<2<1,-1<yS1}U{(@,1):0<2<10<t<1}u 
{(@,-1,8):-1<2<0,0<#<1}, 


Ds = {(,y,):-1<y<10<t<1}U{(1,9,):-1<y<10st< 1}. 


Si riconosce poi (cfr [Ma], paragrafo 5) che 
0,Q = ZoU L2U La. 


La figura 1 mostra la classificazione di Fichera del mantello laterale di Q. La 
base inferiore di Q è costituita da punti di Z2, quella superiore da punti di Z1. 


3. Problema di Dirichlet discreto. 


Determiniamo una griglia G di punti di R° e definiamo su tale griglia un opera- 
tore Lg costituito da differenze finite, che sia consistente con l’operatore differen- 
ziale L. 
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Z1 


Fig. 1 


Osserviamo innanzitutto che la parte del primo ordine nell’equazione di Kol- 
mogorov può essere scritta nella forma seguente 


(3.2) U; — Tuy = - 


dove v è il vettore (0, —z, 1). Utilizzando la formula di Taylor si ottiene 


dun Yy t) _ule + h,y,t) _ 2u(£,y, t) +u(r sa h,y,t) 
(3.3) 0x2 h? 
i h? (du du j 
+3 (Fe +0h,y,t) — pile — 0 e) 
e 


du _ u(r,y—6r,t+6)-u(z,y,t) 602u 

4) Toy) = TLTISTTO Duane , de 
per opportuni 9, 0’, 0” €]0, 1[, e per h, 6 positivi arbitrari. 

È del tutto naturale definire l’operatore Lg come segue: 


Fa (cy -0"65,t+0"8), 


Lau(c,9,)= End uertint—a, 
(3.5) u(r — hy+6z,t— 6) —2u(cr,y+6r,t—6)+u(r+h,y+6x,t-6) 
u(r — hry + 60,6 — 6) — 2u(2,y + 62,6 6) + u(c+h,y+62,t-0) 


h? : 


in modo tale che sia consistente con l’operatore L, secondo la (1.7): 


Leu = Lu+ (h° +6)O(h,6). 
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In analogia con quello che si fa per le equazioni paraboliche, definiamo G in 
modo che abbia passo costante A, lungo l’asse x e A; lungo l’asse #. La definizione 
di Lg data in (3.5) impone allora di scegliere il passo Ay lungo l’asse y nel modo 
seguente: 


Ay = AxAi, 


come avevamo dichiarato nell’introduzione. Poniamo quindi 
G={(jA,FAy,nA;) eR° :j,kne Z}. 


La figura 2 mostra i punti di G che vengono messi in relazione dall’operatore 
La. 


l è i 
e 


Fig. 2 


Nel seguito, per abbreviare le notazioni, scriveremo 


ut = u(jAz,kAy, At), 


1 —1 i -1 
uf lia4i T 20; + Ubi _ VET Viti 
(Ax)? Ai 
Scegliamo per comodità Az = 1/J A; = 1/N, con J,N € N, Ay = 1/(JIN), 
. poniamo 


Loui.k = 


H=GNQ 
Ho ={(jAx,KAy;nA:) eH:n>0,-J<j<J-JIN<k+j<JN} 
0,H = H\Ho 


e consideriamo il problema discretizzato (1.10) 


Lew =0 in Ho 
w= in0d-H. 
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Notiamo che l'equazione LG mette in relazione punti di H solo quando è valutata 
in Ho. Inoltre si riconosce direttamente che 


(3.6) HNdQE9H, dist(P,4,0)<A, VPeò,H. 


4. Stabilità e convergenza. 


In questo paragrafo forniremo delle stime dell’errore tra la soluzione del problema 
(1.5) e la soluzione del problema discreto (1.10). 


Teorema 4.1 [Principio di massimo]. Supponiamo che la griglia G verifichi la 
condizione (1.11) 
Az 
< 
(Ax)? — 


Nim 


e che risulti 


Lew <0 inH, 
(4.1) | ni ; 


w>0 in0,H. 
Allora w > 0 in H. 


Corollario 4.2 [Stabilità]. Supponiamo che valga (1.11) e che w sia soluzione di 


Lcw=0 in Ho 
w=% in0ò,H. 


Allora maxy |w| = maxa,4 |gl. 


Dimostrazione del Teorema 4.1. Procediamo per induzione su n. Per n = 0 il 
risultato è ovvio, in quanto i punti di (Ax, kAy;0) € H sono in particolare punti 
di 0,H. 

Supponiamo ora l’affermazione vera per n. Allora, se (jAz, kKAy, (n+1)A;) € Ho, 
dalla (4.1) segue 


= A Ai n A 
wr > (1-27) utt Wiais +77 ra Wi; 20 
(Ax) (Ax) (Ax) 


per ipotesi induttiva e per la (1.11). Se invece (Az, kAy,(n+1)A:) e 4-H, 
l’affermazione è ovvia. 
Teorema 4.3 [Convergenza]. Sia u la soluzione di (1.5 ), uc la soluzione di (1.10). 


Supponiamo inoltre che Su, du € C(Q). Allora, per ogni e > 0 esiste una griglia 


G, verificante la condizione di stabilità (1.11) tale che 


a TÈ: 
max |u(P) — uc(P)| <e 


Dimostrazione. Fissiamo e > 0 e poniamo 


E =u-ug+e(2- 2°), Et=uc-u+e(2-z?). 
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Osserviamo che, essendo continua e Yu limitata, per la (3.6) possiamo scegliere 
A; in modo che 


|e(P)- uP)|<e 
per ogni P € 9-H. Di conseguenza, nei punti di 0-H risulta 


E > e(2-x?)-e>0, Et>e(2-2?)-e>0, 
mentre nei punti di H si ha 
LeE = Lcu— 2, LeE* =-Lgu-2e. 


Se indichiamo 


dalla (3.3), dalla (3.4) e dalla condizione di stabilità (1.11) segue immediatamente 


2 
[Lsu] < M (w) SET 
in ogni punto di Ho. 
Quindi, se 
8e (Ax)? 
< = 

Ar Si M(u)' Ai 2 ’ Ag AxAi, 
si ha 

LcE  <0 inHo LcE* <0 inHo 

E >0 in0,H, E*>0 in0,H, 


da cui segue, per il Teorema 4.1, 


max lu(P) — uc(P) Se. 
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Soluzione approssimata del problema 


La in Ho: 
con @(x,y,t) 


d_H, 


r 


u=@ in 


per t costante. 


» calcolata 


-3 
28 


t 


A 
(Ax) 


sin(7xx)sin(7y) e 


0.4167 


- 


uzione approssima 


= 

= 
== 

\l 


= 
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ZZZ 
AT DINI 
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5 
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o nu A_ 
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MU 
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PSE 


- 
w 
[=] (°) 
w 
- 
- 
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t= 0.625 


Ù 


\ 


t= 1.042 


